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Abstract 

We study the colour connection structure of R-parity violating decays and pro- 
duction cross sections, and construct a Monte Carlo simulation of these processes 
including colour coherence effects. We then present some results from the imple- 
mentation of these processes in the HERWIG Monte Carlo event generator. We 
include the matrix elements for the two-body sfermion and three-body gaugino 
and gluino decays as well as the two-to-two resonant hard production processes in 
hadron-hadron collisions. 
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1 Introduction 



In the past few years there has been a large amount of interest in R-parity violating (R p ) 
supersymmetric (SUSY) models, motivated by the possible explanations of various exper- 
imental discrepancies, e.g. jn-ljJ. It has become clear that if we are to explore all possible 
channels for the discovery of supersymmetry then R p models must be investigated. For 
a recent review on R-parity violation see JIT^ . 



In the Minimal Supersymmetric Standard Model (MSSM) a discrete multiplicative 
symmetry, R-parity (R p ) is imposed, 



R p = (_l)3B+L+25 ; (1) 

where B is the baryon number, L the lepton number and S the spin of the particle. All 
the Standard Model particles have R p = +1 and their super-partners have R p = —1. The 
conservation of R-parity forbids the terms in the superpotential which violate baryon or 
lepton number [] 

W„ p = l -\ l]k e ab U a LlE k + K, k e ab L\Q{D k + l^ jk e c ^lf c X 2 K + ^H,, (2) 

where i,j,k = 1,2,3 are the generation indices, a,b = 1,2 are the SU{2) L indices and 
ci,2,3 = 1,2,3 are the SU(3)c indices. U (Q l ) are the lepton (quark) SU{2) doublet 
superfields, (D\Tr) are the electron (down and up quark) SU(2) singlet superfields, 
and H n , n = 1,2, are the Higgs superfields. We shall neglect the last term in Eqn.|2| 
which mixes the lepton and Higgs SU(2) doublet superfields. For a recent summary of 
the bounds on the couplings in Eqn. [2] see |TT| . 

This superpotential gives interactions which violate either lepton or baryon number. 
For example the first term will give interactions of a slepton and two leptons which 
violates lepton number. The third term gives an interaction of two quarks and squark, 
which violates baryon number. When combined with the MSSM superpotential there are 
also terms involving the interactions of three sleptons/squarks and a Higgs which violate 
either lepton or baryon number. 

R p is imposed in the MSSM to avoid the simultaneous presence of the second two 
terms in Eqn.|2]. These lead to fast proton decay, in disagreement with the experimental 
lower bounds on the proton lifetime. However in order to guarantee proton stability it is 
sufficient to forbid only one set of these terms. This is achieved for example by lepton 
parity 

(L\ E*) ~{L\&), (3) 
(Q\U\D\H U H 2 ) - (Q^U^D 1 ,^,^), (4) 

which allows the third term in Eqn. ^ but forbids the remaining terms. Thus baryon 
number is violated but lepton number is conserved and the proton is stable. Similarly 
there are symmetries such that lepton number is violated and baryon number is conserved. 

1 It should be noted that some authors choose to define this superpotential without the factors of one 
half in the LLE and UDD terms. This will lead to differences in the Feynman rules, but the results with 
this second convention can always be obtained by taking A or A" to be twice their value in our convention. 
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This also prevents proton decay. Both cases lead to very different phenomenology from 
the MSSM. 

In the MSSM the conservation of R p implies that 

1. Sparticles are produced in pairs. 

2. The lightest supersymmetric particle (LSP) is stable. 

3. Cosmological bounds on electric- or colour- charged stable relics imply that a stable 



LSP must be a neutral colour singlet [12] . 
However, in the case of R p models we can have 

1. Single sparticle production. 

2. The LSP can decay. As the LSP is unstable it does not have to be a neutral colour 
singlet. It can be any supersymmetric particle. 

3. Lepton or baryon number is violated. 

In the MSSM, as the LSP is stable, the experimental signatures of SUSY processes 
typically involve missing transverse energy in collider experiments. However, if R-parity is 
violated, and the LSP decays in the detector, the missing energy signatures of the MSSM 
no longer apply or are severely diluted. It therefore requires a different experimental search 
strategy. In particular, in the B case, where the final state is predominantly hadronic, it 
may be hard to extract a signal over the QCD background in hadron colliders. 

Despite the interest in R p and the potential experimental problems, there have been 
few experimental studies at hadron colliders. The first systematic study of R p signatures 



at hadron colliders was presented in [13 1. More recent overviews of the search potential 



at the LHC and Run II of the Tevatron have been presented in [14, 15|. These studies 



have been limited by the fact that few simulations have been available. In hadron-hadron 
collisions the only available Monte Carlo event generator is ISAJET |16| where the R p 
decays can be implemented using the FORCE command, i.e. the decay mode of a given 
particle, e.g. the LSP, can be specified by hand. However there has been no simulation 
which includes all the decay modes and the single sparticle production processes. 

Here we present the calculations required to produce a Monte Carlo event generator for 
the two-body sfermion and three-body gaugino and gluino R p decay modes as well as all 
two-to-two R p resonant production processes in hadron-hadron collisions. We have only 
included those production processes where a resonance is possible, so for example processes 
which can only occur via a t-channel diagram are not included. However where a process 
can occur via a resonance all the diagrams including non-resonant s-channel and t-channel 
diagrams have been included. We also discuss colour coherence effects via the angular 
ordering procedure (which we describe in detail below), and some preliminary results from 
the implementation of these processes in the HERWIG Monte Carlo event generator JT7 



Details of the implementation of supersymmetric processes with and without R p can be 
found in 0. 

After a general discussion of the angular ordering procedure in the Standard Model 
in Section || we discuss the extension to the R p decays and hard production processes in 
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Section [3]. In Section |] we describe the hadronization procedure which we adopt for the 
fL p processes. We then present some preliminary results of the Monte Carlo simulation 
in Section [|. We have made new calculations of all the necessary matrix elements, and 
include them as an appendix. 



2 Monte Carlo Simulations 

In general a Monte Carlo event generator, for a process involving at least one hadron, 
consists of three parts. 

1. A hard scattering process, either of the incoming fundamental particles in lepton 
collisions or of a parton extracted from a hadron in hadron initiated processes. 

2. A parton-shower phase where the partons coming into or leaving the hard process 
are evolved according to perturbative QCD. 

3. A hadronization model in which the partons left at the end of the parton-shower 
phase are formed into the hadrons which are observed. For processes with hadrons in 
the initial state after the removal of the partons in the hard process we are left with 
a hadron remnant. This remnant is also formed into hadrons by the hadronization 
model. 

We now discuss these three stages in turn. 



2.1 Hard Scattering 

The hard scattering process is described by a matrix element calculated to a fixed order 
in perturbation theory, usually only leading order. The momenta of the particles involved 
in the collision can then be generated according to the matrix element. The Monte Carlo 
event generator then needs to take the results of this perturbative calculation, at a high 
scale, and generate the hadrons which are observed. 



2.2 Parton Shower 

In a scattering process the incoming or outgoing partons can emit QCD radiation, e.g. 
q — > gq and g — > gg, or split into quark-antiquark pairs, g —>■ qq. A full perturbative 
treatment of this part of an event is not possible. (If it were possible it would be included 
in the hard scattering matrix element.) We must therefore make an approximation and 
focus on the dominant contributions in the showering process. The emission of QCD 
radiation is enhanced for (a) Collinear Emission and for (b) Soft Emission. We discuss 
these two below in more detail. Our approximation will consist in focusing on these 
enhanced regions of radiation. We then discuss this in an explicit example. 

(a) Collinear Emission 

If we consider the emission of QCD radiation in the collinear limit then, after az- 
imuthal averaging, the cross section obeys a factorization theorem ||19||. This can 
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be understood as follows, the cross section for a process in which one parton pair 
is much more collinear than any other pair can be written as the convolution of a 
universal, i.e. process independent, splitting function and the cross section for the 
same process where the collinear pair is replaced by a single parton of the corre- 
sponding flavour. Due to this functional form we can then apply this procedure to 
the next most collinear pair in the final state, and so on. We thus have an iterative 
rule which leads to a description of multi-parton final states as a Markov chain Jl9| . 



This can be viewed as an evolution in some energy-like scale, such as the virtuality, 
where a parton at high scale is evolved by successive branchings to a lower scale. 
However, the collinear factorization does not specify what the evolution variable 
should be, i.e. it has the same form for any variable proportional to the virtuality, 
e.g. the transverse momentum. This iterative procedure then correctly resums the 
leading collinear singularities to all orders in perturbation theory ||19|| . 

(b) Soft Emission 

For the emission of QCD radiation in the soft limit, a factorization theorem exists 
for the amplitude of the process. The amplitude for a process in which one gluon is 
much softer than the other energy scales in the process can be written as a product of 
a universal eikonal current and the amplitude for the same process without the soft 
gluon. After we square the amplitude and sum over the spins of the external partons, 
we obtain a result which depends on the momenta of all the external partons. 
It therefore seems unlikely that a Markov description based on sequential parton 
splittings can be recovered. The surprising result PD|,pi] is that, after azimuthal 



averaging, these effects can be incorporated into a collinear algorithm by using the 
correct choice for the evolution scale, namely the opening angle. 




Figure 1: Feynman diagram and colour flow for e + e — > qqg. 



2.2.1 Example: e + e — > qqgg 

We can illustrate this with a simple example, i.e. the process e + e~ — > qqgi, shown in 
Fig. [I]. The semi-classical eikonal current can be used to study the emission of an extra 
soft gluon in this process, i.e. the process e + e~ — > qqgi g2 where the second gluon is much 



4 



softer than the other partons. The matrix element including the emission of the extra 
soft gluon is given by 

M(ki,k 2 ,pi,p 2 ,P3;q) = g s m(k ll k 2 ,p 1 ,p 2 ,p 3 ) • J(g) (5) 

where 

• m(ki, k 2 ,pi,p 2 ,p 3 ) is the tree-level amplitude for the underlying process, 
e + (k 1 )e-(k 2 ) -> q(pi)q(p2)gi(pa)- 

• M(A;i, k 2 ,pi,P2,P3', q) is the matrix element for the process 
e + (ki)e~ {k 2 ) — > q(pi)(l( y p 2 )gi(p3)g 2 (q) , i-e. including the emission of an extra soft 
gluon, g 2 , with momentum q. 

• 3(q) is the non-Abelian semi-classical current for the emission of the soft gluon with 
momentum q, from the hard partons. 

• g s is the strong coupling constant. 

Explicitly in our example, the current, J{q), is given by |21 



J(q) = J b ' M (<7)^, s , (6) 



s=l,2 



where e^ s is the polarization vector of the gluon and 

_ +a ( Pi \ j_ + a +b,q ( P2 \ , • raa'b+a' ( P3 



where t b,q and f abc are the SU(3) colour generators in the fundamental representation 
and adjoint representations respectively. 

Radiation Functions 



We now define the radiation functions of p2[ . This will allow us to express the square of 
the current, Eqn. 0, in a useful way. First we define the dipole radiation function Wij(£l g ), 
which describes the radiation of a soft gluon, with momentum q, from a pair of partons % 
and j. Note that Wij(fl q ) only depends on the direction of q, fl q , and not its energy. 



2 „MO\- f Pi Pj \ _ 2 ( 6 



UJ 2 ^ q ' \Pi-q pyq) ^ \Uj n% HQ 
where 

• u is the energy of the soft gluon, 

• s<; = i$h = 1 - v%Vj cos On, 

1 J 

• £i = l-Vi cos 9i 



(8) 
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• 7» = Ei/rrii = l/y/l - v\ 

• Vi is the velocity of parton i, 

• 9i is the angle between the direction of motion of the soft gluon and the parton i, 

• 9ij is the angle between the partons i and j. 

We can now use the dipole radiation functions to express the current squared for a given 
process in the following form, 



J 2 (q) = -f W(^) (9) 



a 

a; 2 

where C m is the colour factor for the tree- level process, and W(Q q ) the radiation pattern, 
given below in terms of the dipole radiation functions. For the process e + e~ — > qqg the 
colour factor C m = CfN c and the radiation pattern is given by 

W qqg (Q q ) = C A [W qg (Q q ) + W qg (Q q )} - ^W qq (Q q ) (10) 

where Cf = %^ an d Ca = N c are the Casimirs of the fundamental and adjoint rep- 
resentations respectively, with an arbitrary number of colours iV c . This corresponds to 
emission of the soft gluon from colour dipoles, i. e. W qg is emission from the dipole formed 
by the quark and the anticolour line of the gluon, W qg emission from the colour line of 
the gluon and the antiquark and W qq emission from the quark and antiquark. Note that 
the qq dipole is negative which is a problem if we wish to use a probabilistic approach to 
treat the soft gluon radiation. 



The dipole radiation function can then be split into two parts as was done in |22|], i.e. 

w ij (si q ) = wi j (n q ) + wi j (n q ) (11) 

where 

^O-i^) (i2) 

The function W-j has the following properties, 

• After averaging over the azimuthal angle of the soft gluon about the parton i it 
corresponds to emission in a cone about the direction of the parton i up to the 



direction of j, [20. 21 



If the parton i is massive then soft radiation in the direction of the parton is reduced, 
i.e. emission within an angle of order 9 ~ rrii/Ei is suppressed, |22[ . 



In the massless limit it contains the collinear singularity as 9i — > 0, |p2| , but has no 
collinear singularity as 9j — > 0. 

While Wij/ui 2 is Lorentz invariant the individual functions W-j/u 2 and W-j/u 2 are 
not. 
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This allows us to rewrite the square of the current, Eqn. [7], in the following form, using 
these radiation functions, as 

w qqg (n q ) = 2C F (w q % + w*) + c A (w 9 gq + w° q ) 

+N- 1 (W« g - W q % + W& - Wl q ) , (13) 

which should be inserted in Eqn. |9|. This gives our main result in this example. The last 
term in Eqn.[13|, and other terms of this type, can be neglected as it is suppressed by 
1/./V,? with respect to the leading order term, as Cp, Ca oc N c , and is also dynamically 
suppressed as it does not contain a collinear singularity in the massless limit (e.g. the 
singularity in the quark direction cancels between the W^ g and W qq terms.) Thus part of 
our approximation for the parton shower will consist of dropping the l/iV| terms. 



Colour Connected Partons 

We can now define the concept of the colour connected parton. Two partons are considered 
to be colour connected if they share the same colour line. The colour flow, in the large N c 
limit, for the process e + e~ — > qqg is shown in Fig. |I| with a dashed and a solid line. Thus 
the q and g are colour connected and the q and g are colour connected, while the q and 
q are not colour connected. Each quark only has one colour connected partner in a given 
Standard Model Feynman diagram and each gluon has two. Colour connected partners 
are defined at each stage of the iterative parton showering procedure. If the final state q 
were to emit another gluon, g 2 , the new final state q would be colour connected to g 2 and 
no longer to g. g and g 2 would then also be colour connected. 



Angular Ordered Emission and Colour Coherence 

We see from Eqn. ^ that neglecting the final term, using the properties of the function 
Wfj, and averaging over the azimuthal angle of the gluon about a parton, the radiation 
can only occur in a cone about the direction of the parton up to the direction of its colour 
partner. This is shown in Fig. [2|. We can draw a cone around parton one with half-angle 
given by the angle between the momenta of partons one and two. The emission from 
parton one within the cone defined by its colour connected partner, parton two, is called 
angular ordered emission. 

The angular ordering procedure is one way of implementing the phenomenon of colour 
coherence. The idea of colour coherence is that if we consider a large angle gluon it 
can only resolve the total colour charge of a pair of smaller angle partons, and not their 
individual charges. It is therefore as if the larger-angle soft gluon was emitted before the 
smaller angle branchings. There have been a number of experimental studies of colour 
coherence effects. In particular the string effect in e + e~ collisions [p3fl , where there is a 
suppression of soft QCD radiation between the two quark jets in three jet events. There 
have also been studies of colour coherence effects between the initial and final states in 
hadron-hadron collisions, p3,B3. It is now firmly established that event generators that 



do not incorporate colour coherence cannot reliably predict the hadronic final state. 

Although we have averaged over the azimuthal angle of the emission of the gluon in 
both the soft and collinear cases, azimuthal effects, e.g. due to spin correlations, can be 



included p6 |, after the full parton shower has been generated. 
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2.2.2 Non-Planar Colour Flows 



We have explained in an example how the cross section for n + 1 partons factorizes in 
both the collinear and soft limits into a universal splitting term and the cross section 
for n partons. Both of these limits can be implemented by using angles as the evolution 
variable in a Markov branching procedure. We start at the hard cross section, normally 
with a two-to-two process. The maximum angle of emission from a parton is set by the 
direction of the colour partner. We then generate some smaller angle parton, e.g. a gluon 
from a quark. Then we repeat the procedure, i.e. the gluon's colour partner is now the 
colour partner of the original quark, and its anticolour partner the quark, and the colour 
partner of the quark is the gluon. One of the partons will radiate with the maximum angle 
given by the direction of the new colour partner and so on until the cut-off is reached, of 
order 1 GeV, below which emission does not occur. This procedure then resums both the 
leading soft and collinear singularities. 

In processes where there is more than one Feynman diagram it is possible for the 
colour flows in the diagrams to be different. This leads to so called 'non-planar' terms 
from the interference terms, where the colour flows do not match. These are not positive 
definite and hence cannot be interpreted in a probabilistic way for implementation in the 
Monte Carlo procedure. They are always suppressed by inverse powers of N c . A procedure 
must be adopted to split up the 'non-planar' parts of the tree-level matrix element to give 
redefined planar terms with positive-definite coefficients that can be used in the Monte 
Carlo procedure. Such a procedure was proposed in [21] and shown to work correctly for 



all QCD processes. However, as shown in |27|, this is inadequate for MSSM processes 
and hence a new procedure was proposed, which we adopt here. In this procedure the 
'non-planar' parts of the matrix element are split up according to 

\Wfun,i = ^P-|M|L (14) 



|M| 



planar 



where \M\f is the matrix element squared for the ith colour flow, \M\^ lanar is the sum of 
the matrix elements squared for the planar colour flows, and |M|| ot is the total matrix 
element squared. This ensures that the terms are positive definite and the new full planar 
terms have the correct pole structure and sum to the correct total cross section. This can 
then be implemented numerically. 




1. Direction of parton 



2. Direction of the 
colour partner 



Figure 2: Emission in angular ordered cones. 



8 



In this section we have explained how by using a Markov branching procedure we can 
resum both the soft and collinear singularities in QCD. 



2.3 Hadronization 

After the parton shower phase it is necessary to adopt some procedure to combine the 
quarks and gluons into the observed hadrons. This is done in the HERWIG event generator 
using the cluster hadronization model [|28||. This model is based on the concept of colour 



preconfmement. This implies that the invariant mass of pairs of colour-connected partons 
has a spectrum that is peaked at low values, a few times the cut-off used in the parton- 
shower, and is universal, i.e. independent of the hard scale and type of the collision, as 
discussed below and shown in Fig.|TTla. 



In the cluster hadronization model [55], after the end of the parton showering process 
we are left with gluons and quarks. The gluons are non-perturbatively split into light 
quark- ant iquark pairs. The final state then consists only of quarks and antiquarks which 
are, in the planar approximation, uniquely paired in colour-anticolour pairs. These pairs 
of colour-connected quarks do not necessarily have the correct invariant mass to form a 
meson. Instead they are formed into colour-singlet meson-like resonance called 'clusters'. 
These clusters then decay in their rest frame to a pair of hadrons (either two mesons or 
a baryon and an antibaryon) with the type of hadron determined by the available phase 
space. In the original model of |2j| these decays were isotropic in the rest frame of the 



cluster, however in the current implementation of the model ||17|| , the hadrons containing 
the quarks from the perturbative stage of the event continue in the same direction (in the 
cluster rest frame) as the original quark. 

It is reasonable to assume that the low mass clusters are superpositions of hadron 
resonances and can be treated in this way [28]. However, a fraction of the clusters have 



higher masses for which this assumption is not valid and these clusters are first split using 



a string-like mechanism [F2SJ into lighter clusters before they are decayed to hadrons. 

A simple extension of this model is used for hadron remnants. For example in a 
collision in which a valence quark from a proton participates in a hard process, the two 
remaining valence quarks are left in the final state. They are paired up as a 'diquark' 
which, in the planar approximation, carries an anticolour index and can be treated like 
an antiquark. The resulting cluster has baryonic quantum numbers and decays into a 
baryon and a meson. 



3 Angular Ordering in p, p 



In Standard Model and MSSM processes apart from complications involving processes 
where there are 'non-planar' terms pTj the angular ordering procedure is relatively straight- 



forward to implement. However in R p SUSY there are additional complications. 

The lepton number violating processes, which come from the first two terms in the 
superpotential, Eqn. [|, have colour flows that are the same as those which occur in the 
MSSM. On the other hand the baryon number violating interactions, which come from 
the third term in Eqn.0, have a very different colour structure involving the totally an- 
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tisymmetric tensor, e clC2C3 . We look first at the colour structure of the various baryon 
number violating decays which we include in the Monte Carlo simulation and then at the 
structure of the various hard scattering processes. 



3.1 Decays 

From the point of view of the colour structure there are three types of baryon number 
violating decays which we include in the Monte Carlo simulation. 

1. Two-body Fi decay of an antisquark to two quarks or a squark to two antiquarks. 

2. Three-body Fi decay of a colourless sparticle, i.e. a neutralino or a chargino, to three 
quarks or antiquarks. 

3. Three-body J3 decay of the gluino to three quarks or antiquarks. 

In general it is possible to consider for example the decay of a neutralino to three 
quarks as either a three-body decay or two sequential two-body decays, of the neutralino 
to an antisquark and a quark and then of the antisquark to two quarks. If either of the 
two sequential two-body decays are kinematically forbidden, i.e. they can only proceed if 
the internal particle in the three-body decay is off-shell, then we consider the decay to be 
three-body, otherwise we treat the decay as two sequential two-body decays. 

The problem is then how to implement the angular ordering procedure for these three 
processes. We shall consider them using the eikonal current with an arbitrary number of 



colours as was done in Section |2.2.1| for the process e + e~ — > qqg. So in these R-parity 
violating processes this means we need to consider the decay of an antisquark to (iV c — 1) 
quarks and of the neutralino, chargino and gluino to N c quarks. We also have to use the 
generalization to N c colours of the antisymmetric tensor, i.e. e Cl '" CNc . 



3.1.1 Squark Decays 

For the decay of an antisquark to (N c — 1) quarks the leading infrared contribution to the 
soft gluon distribution has the following factorized form. 

M(p ,pi,p 2 , ...,p Nc _i,q) = g s m(p ,pi,p 2 , ■ ■ ■ ,Pn c -i) • J(?) (15) 

where 

• m(po,Pi)P2) • • • ,Pn c -i) is the tree-level matrix element for an antisquark, with mo- 
mentum p , to decay to N c — 1 quarks, with momentum pi, . . . ,pn c -i- 

• M(po)Pi)P2) ■ • • ,Pn c -i', q) is the tree- level matrix element for the decay of an anti- 
squark to N c — 1 quarks including the emission of an extra soft gluon, with momen- 
tum q. 

• J(q) is the non-Abelian semi-classical current for the emission of the soft gluon with 
momentum q from the hard partons. 
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Co is the colour of the decaying antisquark and ci, . . . cn c -i are the colours of the 
quarks. 



Again the current, 3(q), is given by, 3(q) = J2 s =i 2 J^^am where here 



J b >"(q) = - 



Po 



N c -1 



q Xn,ci,...,cjv c 



Po ■ q) c °' c ° 



" + E 



i=l 



Pi 



Pi ■ q 



+6,qj c ,...,<,..., CJ v c 

c,:,c'. 



(16) 



6 and /i are the colour and Lorentz indices of the emitted gluon; t 6 '^ , t b ' qi are the colour 
matrices of the antisquark and quarks, respectively. 

We obtain the soft gluon distribution simply by squaring the current 



3 2 {q) = -C F N c {N c -2)\ 



N c -1 

E 



Pi 



Po-q Pi-q 



N C -1N C -1 

+ EE 

i=l j>i 



Pi 



Po 



Pi-q Pj-q 



• (17) 



This can be expressed in terms of the radiation functions as in Eqn. [|, where here the 
tree-level colour factor is C m = e c o,--; c N c -i e co,---,CN c ~i _ jyj^ -where we have not averaged 
over the initial colours, and the radiation pattern is given by 



W(Q q 



'Nc-1 

E 

i=i 



Po 



Pi 



Po • q Pi-q 



2 7V c -2Af c -l 

+ E E 

i=l j>i 



Pi 



Pi 



Pi-q Pj-q 



We can then re-express this result in terms of the functions given in [22 



2CV 



N,-l N c -1 



(Nc-l) j 



(19) 



i=0 j^i 



This is exactly the same result obtained in [p9], in the context of baryon number violation 
in the Standard Model, except that the massless radiation functions of their paper are 
now replaced by the massive functions here. 

This leads to the following approach for treating the soft gluon radiation from this 
process. The quarks from the decay are randomly colour connected to either the decaying 
antisquark or the other quark. This then correctly treats the soft gluon radiation from 
the decay products. 

In general, the QCD radiation from sparticles, which are here in the initial state, is 
neglected in HERWIG. We would expect this approximation to be valid for two reasons, 
firstly the sparticles will usually have a short lifetime and secondly, due to their heavy 
masses the QCD radiation will also be suppressed unless they have momenta much greater 
than their masses. However for the decays we consider, we can include the effects of 
radiation from the decaying sparticles. This is done by treating the radiation in the rest 
frame of the decaying squark where there is no radiation from the decaying sparticle, 
which HERWIG would not generate anyway. However as stated in Section [2.2. 1| while the 
radiation from individual partons, i.e. W^, is not Lorentz invariant the dipole radiation 
functions are. Hence the total radiation pattern is Lorentz invariant and therefore by 
treating the decay in the rest frame of the decaying particle we correctly include the QCD 
radiation from the decaying particle when we boost back to the lab frame. 
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3.1.2 Chargino and Neutralino Decays 

The charginos decay via the process shown in Fig. ^ and the neutralinos via the process in 
Fig.[|. If we consider the QCD radiation from the decay of a colour neutral object which 
decays, for an arbitrary number of colours N c , to N c quarks, then we see that there is only 
one possible colour flow for this process. The squarks appearing in these processes, q ia , 
can be either of the states a = 1, 2 resulting from the mixing of qn and q^, as discussed 
in more detail in the appendix. 




Figure 3: UDD decays of the x + ■ 

In fact, the colour structure of this process is very similar to that of the squark decay 
and the matrix element in the soft limit can be written in the same factorized form as 
before. Again, we can express the current as in Eqn. ^ where the tree-level colour factor 
C m = e c O'---. c iv c -i e cov,Civ c -i _ jyj anc [ t;h e radiation function is given by 

or n c N c 

w <m = w^t ) Y.TK (2°) 

This result can be interpreted as saying that a quark in the final state should be randomly 
connected to any of the other quarks from the neutralino or chargino decays. 




Figure 4: UDD decays of the x°- 
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3.1.3 Gluino Decays 

The colour structure of the gluino decay is very different from that of the colourless objects 
or the squarks which we have already considered, the diagrams for this process are shown 
in Fig. ^. Again if we consider an arbitrary number of colours, N c , the gluino will decay 
to N c quarks. In this case there will be iV c possible colour flows, corresponding to the 
Feynman diagrams and colour flows shown in Fig. || These different colour flows will lead 
to 'non-planar' terms which must be dealt with. 



U; 



~ \ 
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~j v 

(lja 




d V 
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Figure 5: UDD decays of the g. 

The leading infrared contribution to the soft gluon distribution can be written in the 
following factorized form. 



M(p ,Pi,P2, • • • ,Pn c ; g) = g s ^2mi(po,pi,P2, . . . ,p Nc ) • J»(g) 



(21) 



i=l 



where 



• mj(po,Pi,P2) ■ ■ ■ jPn c ) is the tree- level matrix element of the three-body gluino decay 
for the zth possible colour flow. 

• M(po)Pi)P2) • • • ,Pn c 'i q) is the tree-level matrix element for the three-body gluino 
decay including the extra emission of a soft gluon of momentum q. 

• Ji(q) is the non-Abelian semi-classical current for the emission of the soft gluon, 
momentum q, from the hard partons for the iih possible colour flow. 

Again the current has the form Ji(q) = \^J^(q)£^ s , where in this case 



s=l,2 



( P%_\ f ba>a t a> r, v ,. , ( £_) t b ^ e Cl - c "' 

\P0-q) \Pi-q J dc\ c / c // 

N c 



.c Nc 



^)t b ,t a ,e c 



'22) 



We can write the matrix element squared for this process as 



|M(p ,Pi,P2,---,Piv c ;g)| 2 = g 2 s ^2\m i (p ,p 1 ,p 2l . . . ,p Nc )\ 2 ■ |J<(g)| 5 



i=i 



N c N c 

+9sYl Yl m^m* ■ J 4 (g) • J*(g) 

i=l j=l,jyti 



(23) 
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a. Feynman Diagram 



b. Colour Flow 



Figure 6: Baryon number violating decay of g. 



The procedure of P7|1 , which was described in Section |2] can be used with the matrix 
elements for this process given in the appendix to deal with the 'non-planar' terms. So 
we now consider the radiation pattern of the planar terms. The current can be written 
as in Eqn.|] where here the tree- level colour factor 

Cm = tl cl e Cl ---<--- c ^t b c , /c . e c o- c "- c ^ = C F N C \. (24) 

We have not averaged over the initial colours, and the radiation function is given by 

N c N c 

w(n g )= c A w^ + 2C F w^ + ^ £ wj k + (c A w% + 2C F w& 

N c 

+i< + iv^fciyE (K-n-H)- ( 25 ) 

This planar piece of the soft radiation pattern gives us the result we would naively expect. 
This pattern can be thought of as saying the zth quark should be connected to the colour 
line of the gluino as in an MSSM process and the anti-colour line of the gluino and the 
remaining quarks should be randomly connected in the same way as for the baryon number 
violating squark decay. This radiation pattern, Eqn.|25], also contains pieces which are of 
order l/N% with respect to the leading order pieces that we will neglect as in Section ||. We 
therefore connect the zth quark to the gluino in the standard MSSM way with probability 

\M\ 2 

given by , where \M\ 

2 }uii,i ls gi ven by Eqn.[TJ|. We can then treat the anti-colour 
line and the remaining quarks as a decaying antisquark. 

3.2 Hard Processes 

In addition to the decays which we have already discussed, there are a number of baryon 
number violating hard subprocesses we include in the simulation. All of the colour struc- 
tures of the hard processes that actually violate baryon number have already been dis- 
cussed as these processes are merely crossed versions of the various decays discussed above. 
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However in addition to these processes there are some hard processes that occur via the 
third term in the superpotential but involve no net baryon number violation, e.g. Fig. ^ 
We will therefore only discuss this type of process which cannot be obtained by crossing 
the previous results. 

3.2.1 Resonant Squark production followed by Decay 

As before, we will consider the process in Fig.|7| for an arbitrary number of colours, N c . 



Figure 7: Resonant squark production followed by ^ decay for an arbitrary number of 
colours N c . 

We can write the matrix element for the emission of an extra soft gluon in the following 
form. 

M(pi, . . . ,Pn c -i -h,..., k Nc _i, q) = g a m{px, ■ ■ ■ ,Pn c -i -h,..., fcjv c -i) • J(<?) (26) 



• m(pi, . . . ,Pn c -i '■ k±, . . . , k Nc -\) is the tree-level matrix element for the (N c — 1) 
quarks to (iV c — 1) quarks scattering. 

• M(pi, . . . ,Pn c -i '■ k%, . . . , fcjv c _i; q) is the tree-level matrix element for the (N c — 1) 
quarks to (N c — 1) quarks scattering with the additional emission of a soft gluon 
with momentum q. 

• J(q) is the non-Abelian semi-classical current for the emission of the soft gluon, 
momentum q, from the hard partons. 

• pi, . . . ,pn c -i are the momenta of the partons in the initial state. 

• ki, . . . , kN c -i are the momenta of the partons in the final state. 
Again the current has the form J(g) = \ . J 6 '^*?)^^, where in this case 



Qi(Pi) 




Qi(fci) 



q Nc -iipN c ~i) 



qiV e -i(&jV e -i) 



where 



8=1,2 



J^(g) 





(27) 
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where b and fi are the colour and Lorentz indices of the emitted gluon, respectively. 

We can now obtain the soft gluon distribution by squaring the current. This can be 
rewritten using Eqn.|^ where here the tree-level colour factor is given by C m = N C \(N C —1)\, 
again we have not averaged over the initial state colours, and the radiation function is 
given by 

or Nc-lNc-l 

+uf^ E E + O • < 28 > 

^ c ' %=\ 1=1 

where the partons i and j are in the initial state and the partons I and m are in the final 
state. 

This radiation pattern gives quite an unusual angular ordering procedure. If we con- 
sider one of the quarks in the initial state, this quark should be randomly connected to 
any of the other quarks in the initial state or to the final state, i.e. the probability of 
connecting to a given quark in the initial state and the final state as a whole is equal. If 
the quark is connected to the final state it must then be connected at random to one of 
the final state quarks. Similarly the final state quarks are connected at random to any of 
the other final state quarks or the initial state, and again quarks connected to the initial 
state are then randomly connected to any of the initial state quarks. 



4 Hadronization in p, p 



As we saw in Section [| it is possible to angular order the baryon number violating decays 
and hard processes. It is then necessary to decide how to hadronize these events using the 



cluster hadronization model |28| for a full simulation of these processes. The procedure 
described in Section also works in the MSSM provided that the lifetime of the coloured 
sparticles does not exceed the hadronization time-scale. However some modifications to 
this model are required for Jl p processes. 

In the Standard Model and MSSM cases the colour partner for the colour coherence 
effects and for the hadronization phase are always the same. In the F3 decays and hard 
processes we see for the first time cases where the colour connection for the angular 
ordering and for the hadronization can be different. This is because while the colour 
connection for the angular ordering procedure is determined by the eikonal current, the 
colour connection for the hadronization phase is defined by the colour flow in the leading 
order diagram. When baryon number is conserved these are identical, however when 
baryon number is violated, there are cases where the two are different. 

First we consider the simplest type of decay, i.e. a neutralino or chargino decaying 
to three quarks. The method described in Section |3] correctly implements the angular 
ordering procedure. After the parton-showering phase (and the splitting of remaining 
gluons to quark anti-quark pairs) we will be left with pairs of colour connected partons 
forming colour singlets as well as three further quarks. An example of this is shown 
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in Fig. ^. These three remaining quarks form a colour singlet with baryonic quantum 
numbers, a baryonic cluster. To handle baryonic clusters HERWIG needs the constituents 
to be labelled as one quark and one diquark rather than three quarks, so we randomly 
pair two of them up as a diquark. In our example in Fig. |8] the three quarks in the middle 
together form a colour singlet which is combined into a baryon. 




meson 



baryon 



Figure 8: The Feynman diagrams and colour flows for the hadronization of a J3 neutralino 
decay. 

This procedure is relatively easy to implement in the case of electroweak gaugino 
decays. However it becomes more difficult in the case of the decay of a squark to 
two quarks. If the colour partner of the decaying squark is a particle which decays via a 
baryon number conserving process then the two quarks and the particle which gets the 
colour of the second decaying particle can be clustered as in the neutralino case, e.g. in 
Fig. [5] the Uj, dj, dk should be formed into a baryonic cluster. 

However, if this second particle decays via J3 as well, then the procedure must be 
different, as shown in Fig. [10]. Here, instead of forming one baryonic cluster, we form 
two mesonic clusters. This is done by pairing the d^ randomly with either the d; or d m 
into a standard colour singlet, the remaining quark and antiquark are also paired into a 
colour singlet. This is not the colour connection for the angular ordering but the colour 
connection for the hadronization phase, which is different in this case and determined by 
the colour flow in the tree-level diagram. 

This leaves the case of the gluino decay for which it is easiest to consider the two 
colour lines separately. The colour line should be treated as normal and the anticolour 
line like a decaying antisquark. So if the anticolour partner of the gluino is a Standard 
Model particle or decays via a baryon number conserving MSSM decay mode we form the 
three quarks into a baryonic cluster. However if the anticolour partner decays via a 
mode we then form two mesonic clusters. 

There is one further type of colour flow to be considered which is the production of a 
resonant squark via F4 which then also decays via J3. The correct hadronization procedure 
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Figure 9: Hadronization with one J3 decay. 



in this case is similar to that adopted for the case of two colour connected J3 decays. We 
randomly connect the final state quarks to the colour partners of either of the initial 
state quarks to form a colour singlet. The remaining final state quark can then be paired 
with the colour partner of the other initial state quark. This gives two colour-singlet 
clusters. Again the colour partner for hadronization is determined by the colour flow in 
the tree-level diagram. 




Figure 10: Hadronization with two J3 decays. 

Using the procedures we have outlined above it is possible to hadronize any of the ^ 
decays or hard processes. There is however one potential problem. The cluster model 
is based on the idea of colour preconfinement. In baryon number violating processes we 
see a very different spectrum for the baryonic clusters formed from the baryon number 



violation to that seen for clusters in the hitherto studied Standard Model events. Fig. |TT 
shows the spectra for both Standard Model and baryon number violating clusters. The 
spectrum for the baryon number violating clusters peaks at a much higher mass than the 
baryon number conserving clusters and has a large tail at high masses. This therefore 
means that before these clusters are decayed to hadrons most of them must be split into 
lighter clusters. The spectrum of the baryon number conserving clusters in these events 
has the same spectrum as in Standard Model events. 



Fig. [TT] contains the mass spectrum of pairs of colour connected partons after the 
parton-shower phase and the non-perturbative splitting of the gluons into quark- ant iquark 
pairs. The baryon number conserving clusters, Fig. |Tl~|a, contains all the clusters in e + e~ 
events at the given centre-of-mass energies, whereas the baryon number violating clusters, 
Fig. |IT| b only contains those clusters which contain the three quarks left after all the other 
quarks are paired into colour-singlets from neutralino decays at the given mass. 
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Figure 11: Distributions of the colour singlet cluster masses. The baryon number con- 
serving clusters come from e + e~ events at the given scale whereas the baryon number 
violating clusters come from decays of neutralinos at the given masses. 



The joining clusters in FigfL^a shows the clusters from e + e~ events with one quark 
from the parton-shower of the quark, and an antiquark from the parton-shower of the 
antiquark. The remnant clusters, Fig.|T2|b come from the cluster in DIS events which 
contains the diquark, formed from two of the valence quarks. 

We would expect the baryon-number violating clusters to be heavier than the standard 
baryon number conserving clusters because 

1. The baryonic cluster is formed from three quarks originating from three different 
jets, as also shown in Fig. [11] of the neutralino decay. In normal e + e~ — > hadrons 
events the clusters joining partons from different jets are heavier than the clusters 



which come entirely from partons from one jet, Fig.|T2|a. 



2. The new cluster contains a diquark and in general the clusters containing diquarks 
in for example the hadron remnant in deep inelastic scattering are heavier than the 
normal quark- antiquark clusters, Fig. |P2|b. 

As these clusters are heavier they will be more sensitive to the fine details of the 
hadronization model. In particular these clusters are sensitive to the maximum cluster 
mass before the clusters are split and the details of this splitting mechanism. It is worth 
noting that the same is true of the joining clusters in e + e~ — > 3 jet events, and it is 
precisely these clusters that contribute to the 'string effect', which is well described by 
HERWIG. 
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Figure 12: Masses of the joining and remnant clusters. The joining clusters are from e + e~ 
events at the given scale. The remnant clusters were generated in e~p events with 30 GeV 
electrons and 820 GeV protons with the given scale as the minimum value of Q. 

5 Preliminary Results 

We have implemented the R-parity violating decays and hard processes into the HERWIG 
Monte Carlo event generator according to the algorithms given in Sections |3| and f|. They 
are available in the latest version 6.1 [I7|. Having taken care to implement colour coher- 
ence effects, it is of immediate interest to see whether they have a significant influence 
on observable final-state distributions. To this end we have studied some jet production 
processes and compared the final-state distributions with those from standard QCD two- 



jet events. It was observed in Ref. |25| that certain variables can be constructed that are 
particularly sensitive to colour coherence effects. In particular these variables are sensitive 
to the presence of colour connections that link the initial and final states. To investigate 
these effects for the different colour connection structures of the R p models, we will study 
these variables for two-jet production via resonant sparticle production in hadron-hadron 



collisions. We essentially follow the details of the analysis of Ref. [25]. 
As examples we study the processes 

• ud — > ud via a resonant stau, and dd — > dd via a resonant tau sneutrino, this 
occurs via the coupling A^n. The Feynman diagrams for this process are shown 
in Fig.[l3|. This process involves lepton- number violating couplings but no baryon 
number violating vertices. 

• Resonant squark production via the coupling X"2i2- This leads to resonant down, 
strange and charm squark production. The resonant diagram for this process is 
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shown in Fig. 0. 




Figure 13: Feynman diagrams for ud — > ud and dd — > dd. 

These couplings were chosen to try and maximize the cross section given the experi- 
mental constraints on the couplings. The coupling A'311 has an upper bound, at the 2a 
level, given by [|TT 

A '3n < 0.11 (t^ttt) ■ ( 29 ) 



^100 GeV, 

While the bounds on other LQD couplings are weaker they involve higher generation 
quarks and hence the cross sections will be suppressed by the parton luminosities. 

Similarly for the resonant squark production the couplings that couple two first- 
generation quarks would in principle give the highest cross sections. However, the limits 
on these couplings are so strict that we used A W 2i2, which is only limited by perturbativ- 



ity [11 



We would expect very different results for the variables that are sensitive to the initial- 
final state colour connections for these two processes. The first process only has such 
colour connection in the relatively suppressed t-channel sneutrino diagram, whereas the 
second process has initial-final state colour connections in the resonant diagram due to 
the random colour structure at the ^ vertex. This is the effect we wish to demonstrate 
explicitly below. 

5.1 Resonant Slepton Production 

For resonant slepton production we consider both the signal and the background which 
were generated using the program described in Ref. ]T7|]. The only cut made was to require 
the presence of at least one jet with of greater than 200 GeV in the event. A parton 
level cut requiring the pt of the two final-state partons to be greater than 150 GeV each, 
was used to reduce the number of events we needed to simulate, however this should 
not affect the results. The signal points were generated using the following SUGRA 
parameters, M = 600 GeV, Mi = 200 GeV, A = OGeV, tan/3 = 10, and sgn/x = +. 
At this SUGRA point the right down squark mass is 728 GeV which corresponds to a 
limit of A ; 3n < 0.80. The third generation lepton masses are all close to 600 GeV, 
M fl = 599 GeV, M? 2 = 617 GeV, Mp T = 610 GeV. These processes have previously been 



considered in [15. 50 



The results in all the graphs correspond to the number of events at Run II of the 
Tevatron, with centre of mass energy of 2 TeV, and integrated luminosity of 2fb~ . 
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As can been seen in Fig. Ill, the results for two different values of the coupling show 
that there is a bump in the di-jet invariant mass distribution, Mjj, from the resonant 
particle production for large values of the coupling. 
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Figure 14: Di-jet invariant mass distribution for A^n = 0.7 and A'311 = 0.8. 

We can now study the events around the bump, 580 GeV < Mjj < 640 GeV, in 
the distribution and plot the variables that are sensitive to angular ordering for these 



events [p5 |. These variables depend on the distribution of a third jet in the events which is 
generated in the simulation by the parton-shower algorithm. The three relevant variables 
are: 773, R and a. They are defined in the following way: (1) If we define the jets in the 
order of their Et, with jet 1 being the hardest jet in the event, then 773 is the pseudo- 
rapidity of the third jet. (2) Defining A77 = 773 — 772 and the difference in polar angles 
A0 = 03 — 02, then the variable R = a/ A77 2 + A0 2 . This is the distance between the 
second and third jets in the (77, 0) space. (3) If we first define AH = sgn(77 2 ) A77, we can 
then consider the polar angle in the (|A0|, AH) space. This is a = tan _1 (Aff/| A0|). 

In the analysis in Ref. [^5j additional cuts were imposed in terms of these new variables, 
which we also implement in our analysis 

1. A pseudo-rapidity cut on the two highest pr jets in the event, [771 J , |t7 2 | < 0.7. 

2. Requiring the two leading jets in E T to be back to back ||0! — 2 | — n\ < 20°. 

3. We require the transverse energy of the third jet, Et3 > 10 GeV to avoid background 
from the underlying event. 

4. For the study of a only, we make the additional cut 1.1 < R < tt to avoid problems 
with overlapping cones in the jet clustering algorithm. 
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We can now study the distributions for the signal, the background and the signal plus 
background for the resonant slepton production with coupling A^n = 0.8. There are 
significant differences between the signal and the background for this process. In the rj 3 
distributions, Fig. [1^, instead of a dip in the background at r] 3 = there is a bump in 
the signal. This dip in the QCD background was observed in [25], and is a feature of the 
initial-final state colour connection. In our study it is present in the background, but not 
the signal. 
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Figure 15: The distribution of events in r] 3 for (a) the resonant slepton production, (b) 
the QCD background, and (c) the combination of the two. 



The distribution of events in R, shown in Fig. IIB, is very similar for both the signal 



and background. In fact in the study of |25] all the event generators, even those which 



do not include angular ordering, gave good agreement with the data for this observable. 



The distribution of events in a, Fig.[17|, also shows a difference between the signal and the 
background, with the signal not showing the dip in the middle. This is again an effect of 
the initial-final state colour connection which is present in the background but not in the 
signal. 

As can be seen in all the distributions apart from the disappearance of the dip at 7/3 = 0, 
once the signal and background are added the effect of the signal is minimal. While there 
are differences between the signal and background it is hard to see how cuts can be applied 
on these variables to improve the extraction of a signal over the QCD background. The 
only major difference which can be cut on is the difference in the distribution of a. We 
consider two approaches to increase the ratio of signal to background S/B . 

1. Accept all the events with at least three jets, provided they pass the cuts described 



above from the analysis of [p5| . 

2. Reject all the two jet events and only accept the events with more than two jets 
provided that \a\ < a cut . We apply a cut of a cut = 0.4 for the these jet events. 
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Figure 16: The distribution of events in R for (a) the resonant slepton production, (b) 
the QCD background, and (c) the combination of the two. 
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Figure 17: The distribution of events in a for (a) the resonant slepton production, (b) 
the QCD background and (c) the combined events. 
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These cuts were chosen to maximize S/ B while not reducing S/y B below five. As can 



been seen in Fig. [18| both of these cuts significantly increase the S/B. This can be seen in 
the effect on the invariant mass distribution with the second cut, Fig. [19]. In the invariant 
mass distribution the signal is now more visible over the background. 

This shows that by using the colour coherence effects we can improve the extraction of 
a signal. We would expect obtaining a large S/B to be important for this process because 
we do not have an accurate prediction for the QCD background. However given the 
limits on this coupling this signal will only be visible above the background at the highest 



couplings currently allowed by low energy experiments. In |15|,|30| it was suggested that 
by using the sidebands to normalize the background that resonant slepton production 
could be probed to much smaller values of the couplings. Indeed the S/ \/B numbers 
in Fig. [18] suggest that without any of our additional cuts the signal is visible at much 
lower couplings. However their results were obtained using the narrow-width limit for the 
production cross section and did not included the effects of QCD radiation. Our results 
suggest that after including these effects the signal will only be visible for large values of 
the coupling. It may be possible to use the sidebands which we have removed with our 

1B0| 



cuts to normalize this background, as in 



to improve the extraction of the signal. 



However this may not be possible due to the increased width of the resonance, Fig. [14], 
due to QCD radiation. The situation will hopefully improve with the availability of a 
next-to-leading order calculation for the QCD background. 
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Figure 18: Effect of the cuts on the angular ordering variables as a function of A'311. 



5.2 Resonant Squark Production 

The cross section at the Tevatron for resonant squark production^ is much lower than the 
resonant slepton production due to the reduced parton luminosity. 

(This will be reversed at the LHC.) It is therefore unlikely that an excess of events 
can be seen over the QCD background. It is however still interesting to look at the 

2 Resonant squark production via UiDjDk has previously been considered in fjlf . 
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Figure 19: Invariant mass distribution for A'sii 
angular ordering variables. 



0.7 and A'311 = 0.8 after cuts on the 



distributions of the angular ordering variables. The signal point was generated using the 

following SUGRA parameters, M = 500 GeV, Mi = 200 GeV, A = OGeV, tan/5 = 10, 

2 

and sgn/x = + and X'212 = 1- At this point the squark mass § g = 601 GeV. As can 
be seen in Fig. |2(| there is now less difference in the shape of the distribution between the 
signal and the background (c.f. Figs. [15], [Tj| [I7p. The resonant squark production shows 
a dip at r/ 3 = and a rise as a — > |, which is due to the colour connection between the 
initial and final states. The effect is slightly less than for the QCD background as there 
are combinatorially fewer such connections. 

The fact that the final state distributions of the resonant slepton and resonant squark 
production processes are so different, despite the identities and kinematics of the jets 
themselves being so similar, clearly shows that colour coherence plays an important role 
in determining the properties of R-parity violating processes. Even if this is not used as a 
tool to enhance the signal, it is likely that it will affect the efficiency of any cuts that are 
applied, so it is essential that any experiments looking for R-parity violating processes 
take into account colour coherence in their simulations of the signal. 

Even if R-parity violating hard processes were added to ISA JET [TjJ , it would not be 
expected to describe the final state well, as it is based on the incoherent parton shower and 
independent fragmentation models. Thus, in our case for example, the resonant slepton 
and resonant squark processes would have very similar properties. It is worth noting that 
ISAJET gives a poor description of the CDF data ]25| on 773, R and a in standard QCD 
two- jet events. 
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Distribution of rj 3 



Distribution of R 



Distribution of a 




Figure 20: Distributions of the signal for resonant squark production. 

6 Conclusion 

We have presented a procedure for implementing colour coherence effects via the angular 
ordering procedure in R parity violating SUSY models. We find that the baryon number 
violating processes have a random colour connection structure for angular ordering. In 
these processes we see for the first time differences in the colour partners for the colour 
coherence effects and those used with the idea of colour preconfinement for hadronization 
in the cluster model. 

A full set of decays and hadron-hadron cross sections have now been implemented 



in the HERWIG Monte Carlo event generator |T^,|T8[. The first preliminary results for 
these processes show that the inclusion of colour coherence is important and that for some 
processes we can can use the colour coherence properties of the processes to help extract 
an R parity violating SUSY signature. 

The availability of a full simulation should allow a more detailed experimental study 
of these processes for the first time. 
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A Conventions 

Here we present all the matrix elements for R-parity violating two- and three-body decays 
of sparticles as well as the the matrix elements for single sparticle production via 2 — > 
2 scattering processes. We disregard those possibilities where the sfermion resonance 
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kinematically is not probed, e.g. 



dj + d k 



Vi 



h + z°. 



(30) 



We also do not consider the processes generated from quark-photon scattering as discussed 
First we present the decay matrix elements of the sfermions, neutralinos, charginos 



in 
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and the gluino. We will then give the matrix elements for the cross sections most of which 
can simply be obtained by crossing the various decay matrix elements. Throughout we 
allow for more than one R p coupling to be non-zero. 

We follow the conventions of [|3^,|33j for the neutralino and chargino mixing matrices 



and the convention of |35[ for the mixing of the sfermions. For the current eigenstates 
q l L R and the mass eigenstates q\ 2 the mixing is 



(31) 



We denote the mixing matrix above as Q l j k where i — u, d, s, c, b, t is the quark flavour in- 




dex. The analogous slepton mixing matrix is denoted Lj k , where i = e~ , u e , fjT, u^, 
is the lepton flavour index. We neglect inter-generational sfermion mixing. As we do not 
consider the right-handed neutrino we also neglect the lepton mixing. We give the formu- 
lae below for general generation indices. However, in HERWIG, the mixing for the first 
two generations of sleptons and squarks is not included as it is expected to be small. 

In order to simplify the notation for the matrix elements we introduce the following 
functions 



R(a, m 



be) 



K c - M?y + r?Mf 



2 ; 



(32) 



S(a, b, m 2 cd7 mj f ) = R(a, m 2 cd )R(b, m 2 ef ) [(m 2 cd - M 2 )(m 2 ef - M?) + r^M.M^ .(33) 

Here m\ c = (pb + p c ) 2 , and M s , T s are the mass and the decay width of the sfermion a, 
respectively. The various terms in the matrix elements can be more easily expressed in 
terms of 



*(5,1,2,3) 
T(5,l,2,3) 

$(5,6,1,2,3) 



R(a, m\ 2 ) {m\ 2 -m\- m 2 2 ) 

[(a 2 (a) + b 2 (d)) (M 2 + m\ - m\ 2 ) + 4a(a)6(a)m 3 M ] 
S(di, a 2 , m\ 2l m\ 2 ) (m\ 2 — m\ — m^) 
[(a(ai)a(a 2 ) + 6(ai)6(a 2 )) (Mg + ml - m 2 12 ) 
+2 (a(5i)6(a 2 ) + a(a 2 )6(a x )) m 3 M ] , 



(34) 



(35) 



S(a,b, m 2 12 ,m 2 23 

+m 1 M b(d)a(b) (m 
+m 3 M a(d)b(b) (m 



mim 3 a{a)a 



(b){ 



in 



12 



m 23 - m 1 



m q 



m 



m, 



m: 



m: 



+b(a)b(b) (ml 2 ml 3 — m\ml — M^m^) 



(36) 



Here a\ and a 2 are the mass eigenstates of the relevant SUSY particle. The functions a 
and b are gaugino-sfermion-fermion coupling constants and are given in the tables below: 
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Table for the neutralino, Table ^ for the chargino and Table [| for the gluino. The cou- 
plings are defined such that a(c*) = 6(c), and b(c*) = a(c). In all the above expressions 
M is the mass of the decaying sparticle and 1, 2, 3, are the decay products. 



B Decays 
B.l Sfermions 

Here we present the matrix elements for the two-body sfermion decays including left/right 
mixing. In general the spin and colour averaged matrix elements have the form 

\M(a - b, c)| 2 = C a bc {Ml -ml- m% (37) 

where Cg c is the colour factor and the coupling for the process. These factors are tabulated 
for the various sfermion decays in Table [I]. N c denotes the number of colours. 



Operator 


Process 


Colour Factor and Coupling C£ c 


LLE 




\\ 1 2 1 r 2 j — 1 1 2 
\ A ijk\ \^Xa \ 


LLE 


&ka > V $j 


IX |2|r2Jfe-l|2 
\ A ijk\ \ ±J 2a \ 


LLE 




1 X I 2 

I ^ijk | 


LQD 


— * Ujd k 


N c \K jk \ 2 \L?- l \ 2 


LQD 


Vi — > djdk 


N c \K 3 k? 


LQD 


djot v%d k 


\\) |2|/o2j-l|2 
\ A ijk\ IVla 1 


LQD 


Uj a ► 6j d k 


|\/ |2|/^2j|2 
l^jjfel IVlal 


LQD 


dka ^ l^idj 


1 y |2| / o2fc-l|2 
I'Vjfcl IV2a 1 


LQD 


dka * Uj 


1 y |2| / o2fc-l|2 
l A yfcl IV2a 1 


UDD 


Ui a > djdfc 


(N c -l)\\K jk ?\Qf a \ 2 


UDD 


dka * Uidj 


(N c -l)\\K jk \ 2 \Qt k - l \ 2 



Table 1: Coefficients for the Scalar Decays. 



In all these terms the Roman indices represent the generation of the particle and the 
Greek indices the mass eigenstate of the sfermions when there is mixing. The decay rate 
can be obtained by integrating over the two body phase space. This gives 

tv h \ |M(a ^ 6,c)| 2 p cm 

Via— > b, c = , 38 

V ' ; 87rM a 2 ' V ; 

where p cm is the final-state momentum in the rest frame of the decaying particle 

PL = 4^ [ M l - K + m c ) 2 ] [M 2 a - (m b - m c ) 2 } . (39) 
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B.2 Neutralinos 



The total three-body decay rate of a photino was first computed in [M in the limit 
where the sfermion is much heavier than the decaying photino and assuming massless 
final states. In |37J the general photino matrix element squared was given, allowing for 
the computation of final state distributions. In |38|,|39|j this was extended to the general 
case of a neutralino. In []39| arbitrary sfermion mixing was included as well. We have 
recalculated the rates with only left/right sfermion mixing (neglecting intergenerational 
sfermion mixing). We use a different convention both for the R p superpotential and the 
MSSM Lagrangian, which is more appropriate to the implementation in HERWIG. The 
LLE, LQD and UDD decay modes are shown in Fig. [21], Fig. |22| and Fig. (| respectively. 
There are four decay modes 



1. x° 

2. x° 

3. X° 

4. x° 



^djd k , 

£+Ujd k , 

Uidjdk, 



as well as their complex conjugates, since the neutralino is a Majorana fermion. 






Figure 21: LLE decays of the x- 
The spin and colour averaged matrix elements are given below. 3 



\M(xl 



\ 2 
A ijk 



3 k 



1 1 ? ^ki 



+ ^\L?- 1 \**(£ ja ,i'iM) 



a=l,2 



kai u ii c ji c k) 



a=l,2 



+2L%- 1 Lli- 1 ?(£ j , i/,, 4, tj) + 2Lg" 1 Lg- 1 T(§ , 

- 2^r^fea, vu 4, £j) - E 2^-^(4 



fej u ii Ljj Lk) 



kai ii ii ^jl ^k) 



a=l,2 



o=l,2 



3 We have a slight disagreement with |3^] concerning the sign of the width of the sfermions. This is 
numerically insignificant since when the sfermion is on-shcll HERWIG treats this as a two-body decay. 
The authors of |39| agree with our signs. We thank Paolo Gondolo for discussion of this point. 
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-EE 2L ia 1 ^- 1 m„ e ja , ij, Vu 4) 
0=1,2/3=1,2 



|M(x?-MA)| 2 = 

*%k N c 4, 4, f») + E \Qia l \ 2 ^(d jw u h d k , dj) 

a=l,2 

+Ei^rT^(4 a ,^,^,4) 

a=l,2 

+2Q 2 r 1 Q?r 1 r(4 ^, 4, dj) + 2Q 2 2 k 1 - 1 Qlt 1 r(d* k , v % , dj, d k ) 

- E 2 ^ir 1$ (^«' ^> U i> 4, 4) - E 2 ^2a _1 *(4a ? ^, 4, 4) 



a=l,2 



a=l,2 



- £ £ 2Q£" 1 Q*" 1 <I>(4/j, d,- a , 4, I/*, 4) 

a=l, 2/3=1,2 



|M(x?-^4)l 2 = 



\' 2 ijk N c 



J2 \ L Ta T*(4*,%,4,4) + ^ iQlU^iUjaJi^^Uj) 
t=l,2 a=l,2 

+ E IG£~T*K«, u v 4) + 2L*" 1 L*- 1 T(4 Uj, 4, 4) 

a=l,2 

+2Q 2 iQ 2j 2 T(^-, 4 4, + 2g^- 1 g^-i T (d*, 4, 4) 

- £ ]T 2L?- 1 g 2 ^(^ /3 , 4 a , 4, 4, uj) 
0=1,2/3=1,2 

" E E ^la'Qff 1 ^, L, 4, «,-, 4) 

0=1,2/3=1,2 



-EE (4/3, %a, Uj, 4, 4) 

0=1,2/3=1,2 



|M(x?-^44)| 2 = 

E 1^201^(^4,4,^) + iQt'^idj^u^d^dj) 

La=l,2 a=l,2 

+ E IQfa'^WL, u h dj, 4) + 2Qf 1 Qf 2 T(u*, dj, 4, Ui) 



0=1,2 



+2QH 1 Qli 1 T(d*, Ut , 4, 4) + 2g^- 1 g^- 1 T(4, Ml , 4, 4) 
-EE 2 Q 2 ^ l Q%~ 1 ^, u* a , Ui , 4, 4) 



a=l,2/3=l,2 
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a=l,2/3=l,2 



a=l,2/3=l,2 



(43) 







Figure 22: LQD decays of the 



When the neutralino mass matrix is diagonalized it is possible to get negative eigen- 
values in which case the physical field is 75% rather than \. This then changes the sign 
of some of the coefficients in Table|2|: the coefficients a(c) change sign, and hence the 
coefficients b(c*) also change sign. 

The partial widths can be obtained from these matrix elements by integrating over 
any two of m\ 2 , JT&23 an d m i3- The partial width is given by ]40| 



r(0-l,2,3) 



(2tt) 3 32M 3 



*? 2 



dm\ 2 



rain 



dmlz\M\\ 



(44) 



where 



12/ max 
2 ) 

121 min 
2 



(M - m 3 ) 
(mi + m 2 ) 



(™2 3 ) 



VMX 



[m: 



23/ 



(E* + Elf - (VEf-ml + y/Ef - ml) , 



E* 2 



[m 



12 



ml + ml) /2m 12 and E* = (M 2 



m 



12 



m|) /2m 12 are the energies of 



particles 2 and 3 in the m\i rest frame. 
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Coefficient 



a{pi) 

m 

a{t ia ) 



2 cos 9 W 



mi 



gN [3 T 2i-l 



« 2M W cos (3 



.2«-l I a AT' _ gsin 2 6 w N'i2 
cos W 



6(4 



i-l T 2i-1 ( cM i i gN' l2 (\-sin 2 e w ) 



mn _9Nl3_T2i-l _ tzi-l I „ at/ i 

H i2M w cos/3 2a ^la I 01 '1 ~ COS tfu 



6(4 



Q2a + Via ee d l\ n — -g- 



* 2M w cos/3 



a(ui 



m 



gNu rP-i n 2 i 



U J 2M W sin 



QZ - QZ ee u N' 



ji ge u sin 2 W N', 



cos 9 U 



b{ui 



m 



gNu r\2i 



u i 2A/^sin/3 



Qfa + Qfa ee u N' n + 



gN'i 2 (\-eu sin 2 6 W ) 



Table 2: Couplings for the Neutralino Decays. 



B.3 Charginos 



Most of the chargino fl p decay rates have already been calculated in || in the case of 
no left/right mixing for the first two operators in the ftp superpotential. We recalculate 
these rates with left/right mixing. First we consider the LLE decays of the chargino. 
There are three possible decay modes: 



1Ai 



2. X H 
3- X A 



The Feynman diagrams for these decays are shown in Fig.|23|. The spin averaged matrix 
elements are given by 



\M(xt ^ Pii+iy k )\ : 



9 2 X 2 



a=l,2 



(45) 



9 2 X 2 



5>£-T*(£ 



lai i ^ki^i 



)T l 2 1 1 ^(i jn .u,.if. 



o=l,2 



a=l,2 
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Figure 23: LLE decays of the \ + ■ 



2-^11 T(4j Vji 4> V% 

) + 2L l { Ly 2 T(£j,Ui,£ k ,Uj) 



+ EE ^L'S'Hijp, 4*, u it 4, vj) 

a=l,2/3=l,2 



(46) 



|M(# 



^[*(j> i) ^,4,^)+*(i?i,44,^; 

+2$(? i ,? i ,44^)] 



(47) 



We go beyond the results of |2j to include the decay \ + 

We now consider the LQD decays of the chargino. There are four possible decay 



modes: 




i. x + - 


— > z/idjUk, 


2. X + - 




3. X + - 




4- X + - 


— > ^iUjd k . 



The Feynman diagrams for these decays are shown in Fig. ^4[ The spin and colour aver- 
aged matrix elements are given below 



|M(^ + -> M^)l 2 = 



o=l,2 



(4£ 



|M(x+^£+%« fe )p 



«=1,2 



(49) 
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Figure 24: LQD decays of the \ + ■ 



|M(x+ -> £+4-4 



^{pi.d^dk.ii) + ^ |Q 2 ^| 2 #(w jQ ,44,4) 



0=1,2 



+2Q^Q^T({t i , £j, 4, d,-) + 2 ^ z> i; 4 4, dj 



(50) 



o=l,2 



s 2 a /2 ^ c 



|M(# -> ^4)| 2 = 

1=1,2 0=1,2 
h2Lff 1 L^Tfe u„ 4, Pi) + 2Q 2 f 1 Q 2 f 1 r (4, ^ 4, Uj) 

1/3 



V V ^ Lf a 1 Q 2 1 J 1 ${d j p, I ia , Vi, 4, Uj) 



o=l,2/3=l,2 



(51) 



We now come to the baryon number violating decays. We do not assume here that 
there is only one non-zero fl p coupling. This means that more than one coupling con- 
tributes to these decays. It may seem that this will only matter in the case where more 
than one A" coupling is taken to be non-zero, however there can be more than one dia- 
gram even with only one coupling non-zero, e.g. A'/ 12 will give two diagrams for each of 
the decay modes. In this case one of these diagrams is obtained from the other simply by 
crossing the identical fermions in the final state. 

There are two possible decay modes: 



2. t 



UiUjdkj 
didjdk. 
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The Feynman diagrams for these decays are shown in Fig. [3]. The spin and colour averaged 
matrix elements for these processes with left /right sfermion mixing are given below. 

|M(x+ -> u iUj d k )\ 2 = 
9 2 N C \ 



2(1 + 5, 



x' 2 jik E I^T*(4. u v 4, ud + x% £ |cS, _1 l a *(^a. «i. 4, %) 



a=l,2 



a=l,2 



-2^"?ifcQ2i 1( ?22 ^{d^u^d^Ui) + 2X" 2 ijk Q 2 2 \ X Q% ^(d*,^, d k ,Uj) 



+2\" ijk \" jik X^L 1 ^^ X ®{d*ji3->d* a -: u i-:dk-,Uj) 

a=l,2/3=l,2 



(52) 



a=l,2 



2(1 + Sij + 5 jk + 5 ik ) 
+X '% Yl \Q 2 2e\ 2 ^(v* ja) di, d k , dj) 

a=l,2 



a=l,2 



+2A'^QMq|T(u*, di, 4, d 3 ) + 2A%.QMT(^ a , dj, d,-, d k ) 

— 2\" ijk \"j ki X X ^a^t^^i/?' Uia> di, d k , dj) 
a=l,2/3=l,2 

-2\"i jk \" kij X X ^aaQ^^C^, di, dy, dfc) 



0=1,2/3=1,2 



-2X"j ki X" kij X X ^aQip^i^lp^ dy, G?j, d fc ) 

a=l,2/3=l,2 



(53) 



The coefficients in the chargino matrix elements are given in Table^ and the partial 
widths can be obtained by integrating the matrix elements in the same way as for the 
neutralino decays. 

Again when the chargino mass matrix is diagonalized negative eigenvalues can be 
obtained and the fields must be rotated. This means here that the coefficients bfa), 
b(uia), and 6(dj a ) change sign if the chargino mass is negative. 

B.4 Gluinos 

These decay rates are calculated here with left/right mixing. There are three possible 
decay modes, two via the LQD operator and one via the UDD operator: 

1. g — ► ^djd k , 

2. g — ► £+Ujd k , 
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Coefficient 




Coefficient 




a{£ ia ) 





b(L) 


j2i-\Tj U l2 I% a 1 m ei 
^lot u ll V2M w cosP 


a{vi) 


\/2M w cos p 


Kh) 


V* 
v ll 


a(u ia ) 


\/2% cos/3 


b(u ia ) 


•\r*r\2i '""i^S'^lL 
v lV°tla y/2M w sin/3 


a(dia) 


V2M W sin/3 


b{d ia ) 


Via "11 ^/2M wCOS /3 



Table 3: Couplings for the Chargino Decays. 



3. g — ► UiUjd 



jU k . 



Since the gluino is a Majorana fermion the charge conjugate decay modes are possible as 
well. The Feynman diagrams for these processes are in Fig.|25| and Fig.|5|, respectively. 
The spin and colour averaged matrix elements with left/right sfermion mixing are given 
below. 



\M(g -> Uidjd k 
A' 2 q 2 

7 ' ^\Q?- 1 \ 2 *(d ja ,v i ,d k ,d j )+ J^lQfa'^WL^i^^d, 

a=l,2 a=l,2 

+2Q 2 i r 1 QH- 1 T(d j , v iy 4, dj) + 2Q%- 1 Q 2 2 t 1 T(dt, v u dj, d k ) 



- / ] 1 Q% l ${d* k p,d ja ,dj,i> h d k ) 

a=l, 2/3=1,2 



(54) 



\M(g -> Quid* 



Q = l,2 



Q = l,2 



^QlotOfitp l ^{d>kpi u jaiUj,ii,dk) 

a=l, 2/3=1,2 



(55) 



|M(sf-> 

\'% k (N c -l)\ 



E l^l 2 *( fi to. ^ 4, ui) + 2g|g^T(«*, d i; 4, ui) 



=1,2 
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+ ^2\Q 2 2a 1 \ 2 ^(d* ja ,u i} d k} d j ) + 2Q 2 2 { 1 Q% l T(d*, u h d k , dj) 
0=1,2 

+ E iQl^fWla, u i: d jt d k ) + 2Q^ 1 Qg- 1 T(4 u h d jt d k ) 
1,2 

;3tE E^^W^Xo^^) 



0=1,2 



a=l,2/3=l,2 



iV, 



^jEE ^K/3, U* a , Ui , dj, d k 



"o=l,2/3=l,2 



+ 



N, 



(56) 

"0=1,2/3=1,2 

The coefficients for these matrix elements are given in Table f|. As the gluino mass is not 
obtained by diagonalising a mass matrix it cannot be negative. 



Coefficient 




Coefficient 




a(u ia ) 




b(u ia ) 


-Qfa 


O'idia) 


Ql- 1 


b(dia) 





Table 4: Couplings for the Gluino Decays. 



1A 



djo X 



U; 




Ujo 




dko 




dko 




Figure 25: LQD decays of the g. 



C Hard Processes 

All of the single neutralino, chargino and gluino production cross sections can be obtained 
by crossing from the decay matrix elements we have already presented in Section This 
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crossing will lead to the invariants m 2 2 , m 2 3 , and m 2 3 being replaced by the usual invariants 
s, t, and u. There is an overall sign change due to exchanging fermions between the initial 
and final states. It should also be remembered that the decay matrix elements given have 
been averaged over the spin and colour of the initial particle. The cross-sections for the 
remaining processes are presented below. In all cases the formulae have been averaged 
over the initial spins and colours. The initial state masses have all been set to zero, except 
where they appear in a coupling constant. In t- and w-channel fermion propagators the 
fermion masses have been neglected as well. 



C.l LQD Processes 

C.l.l Resonant Slepton followed by Weak Decay 

There are three processes which can occur via the production of a resonant slepton followed 
by a weak decay of this slepton. These are: 

1. djd k — l?w-, 

2. Ujd k — £*W+, 

3. ujd k — ► f*Z°. 

It should be noted that we have not included processes where the resonance is not acces- 
sible, e.g. Ujdk — ► t~2 Z°. 



\M{djd k - £* a W-)\ 2 = 



s(s-M?.)R(u it s) 



2u 



(m 2 w (2M\ -u)+u(§-Mf )) 

\ HQ HQ / 



(57) 



M(ujd k -> i/TW 



n 2 \' 2 
9 A ijk 



+^ (2M^(ui-M?M^)+u 2 s) 

{M^Ml - u) + u(s - M?J) 



,a=l,2 



|L^-1|2 S( $ _ M 2 j fl(< . ajS) 



2u 



(5S 



\M( Uj d k -> t n Z 



* r 7 0\\2 



N c M'z cos 2 e w 



E l^r 1 ! 2 !^ 1 ! 2 ^^, s 2 ) + 2Lfr 1 Lfr 1 zgz%?i? em s(e il ,i a , s, s) 



,a=l,2 



(2M\{ui - M\ M 2 Z ) + w 2 s) 



+ 



+ E 



a=l,2 



(2M 2 z {ut - M\ Ml) + t 2 s) 
M 2 (2M| i -«)+^( S -M| i )) 
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Q = l,2 



(Ml(2M| i -t)+£( S -M| i )) 



2\Lf- 1 \ 2 Z. 



ut 



^af2M|(M| -t)(M? -u)-stu) 



where in all the above equations 

P 2 c™ = [* - ( m i + m 2) 2 ] [s - (mi - m 2 ) 2 ] , 
and mi, m 2 are the masses of the final state particles. 



(59) 



(60) 



Slepton Couplings 






zf 


I(L^ 1 L^ 1 -2sin 2 M Q/3 ) 


Squark Couplings 




\{-Q¥ a Q% + 2e u sin 2 6 w 5 a p) 


zf 




Quark Couplings 


Z u L 


-\ (l-2e u sin 2 ^) 


Z d L 




Z u R 


2&u sin 6 W 


Z d R 


\e d sin 2 6 W 



Table 5: Couplings of Sleptons, Squarks and Quarks to the Z°. 



C.1.2 Resonant Slepton followed by R p Decay 

There are four process it which we can produce a resonant slepton via R p which then 
decays back to Standard Model particles via a R p decay. These are: 



1. 


djd k — 


-> did m , 


2. 


Ujd k — 


-> uid m , 


3. 


djdk — 


~* A An) 


4. 


Ujdk — 





The first two of these process only require non-zero LQD couplings whereas the second 
two require both non-zero LQD and LLE couplings. The matrix elements are presented 
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below for an arbitrary number of non-zero R p couplings 

\~M(djdk — > dld m )\ 2 = — X'ijkX'ilm^'njk^'nlmS(l'i,l'n, S, s)§ (§ 



m A . - m„ 



_i _ x' x' x' x' 

~r . 7 /v jj' 7 * ikm^ njl" nkm 



\M.(ujdk — > ttzc^ 



(m| ~t)(mg m -£) 
(i-M?)(i-M?J 

' X' X' X' X' I 7- 2i— 1 1 2 1 r2n-l|2 

. / y / A ijk A ilm A njk A nlml^la \ \ Lj \fi \ 



(61) 



i,n=l,3 



a,/3=l,2i,n=l,3 



J 



(62) 



\M(d j d k ^£ l C)\ 2 = X 'ijk X 'njk\lmKlmS(i>i,i>n,S,s)s (s ~ 



2 2 \ 



i,ra=l,3 



(63) 



\M.{ujdk — > 



1 V X' X' X X I j2i-l\2\ j2n-l\ 

^ 2-^ A ijk A njk A ilm A nlm\-^l a \ I 

' '' r a,/3=l,2 i,n=l,3 

S(e ia ,e n p,s,s)s (s -m|J 



(64) 



C.1.3 Resonant Slepton Production followed by Higgs Decay 



There are a number of processes which can occur via the production of a resonant slepton 
which can then decay to either a neutral or charged Higgs: 



1. ujdk 



^ho/Ho/Ao, 



2. djd k — 4H-, 

3. Ujd k — ► z>*H+. 

As we only include processes where there is a possibility of a resonant production mecha- 
nism, the process djd k — > z/fh /H /A is not included. For the same reason we also have 
not included the processes Ujd k — ► £* L h /H /A for the first two slepton generations. This 
is because HERWIG does not include mixing for the first two generation sleptons and the 
initial state only couples to the left-handed slepton. The process Ujd k — ► £* 2 ho/Ho/Ao is 
also not included for the third generation (i = 3) as there is no accessible resonance. 

The matrix elements for these processes are given below. Since the matrix elements 
have the same form for all the neutral Higgs processes we use the notation Hj, where 
/=1,2,3 is h , H and A . 



\M(djd k 



,2 V 2 . 



-M2 



ff 2 V 



4N C 



\H c ~~ e 

1 Vli, 



2 sR(ui, s) + 



4\Lf 



2i ~ 1 | 2 |Q|| 2 



ut 



M L M »-) 



(65) 
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m Ui cos a 




rrid- s i n a 




2M W sin P 


2M W cos P 


u! 


m u . sin a 


A 2 


Trid i cos a 




2M W sin /3 


2M W cos P 


uf 


ra Ui cot P 
2M W 


Df 


m di tan 
2M W 




m Ui cot P 
2V2M W 


D i 


vn,d i tan /3 



Table 6: Higgs Couplings to quarks. 



im(^4^;# + )P 



4N C 



E i^ _1 i 2 itf&j 2 s*&.) +2L-- 1 ^ 2 r 1 ^ i ^ 2 ^ 1 ,4,s,s) 



.a=l,2 



(66) 



|M( Mj 4-^o)! 2 



9 2 A' 2 



47V C 



,a=l,2 
4 



4/3 #0 



(67) 



The couplings involved in the various processes can be found in Tables^ and [7|. 
C.2 UDD Processes 

C.2.1 Resonant Squark followed by Weak Decay 

There are four process which can occur via the production of a resonant squark followed 
by a weak decay of this squark: 



1. djd k 

2. djd k 

3. Uidj - 

4. uidj - 



d*W- 



u£W H 



b*Z 



Again we do not include processes where the resonance is not accessible, i.e. djdk 
and Uidj — > b^Z . So the matrix elements for these processes are given by 



t*Z° 



\M(djd k ^d* W-)\ 2 
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Coefficient 



H 



1 



M z sin(g+ / 3) r j 2i-l j 2i-l / 1 
cos6U L-^la ^1/3 1 2 



sin 2 ^)+ sin 2 ^L^- 1 ^- 1 ] 



m; . sm a 



r r 2i-l r 2i-l _| r2i-lr2i-l] 
L^la ^1/3 ~+~ -^2a ^2/3 J 



Mw cos /3 

^ (^cosa + A ei sin a) [L^L*" 1 + L^Lfr 1 ] 



2M W cos /3 



H l In 



M z cos(a+/3) r 7 -2i-lr2i-l/l ■ 2g \ , - 2 r 2i-l r 2i-l] 

m ej cosa \j2i-lj2i-l , r 2i-l r 2i-l] 
M^cos/3 l^la ^1/3 ^ ^2a ^2/3 J 

(jjl sin a - A ei cos a) [ifc 1 !^ 1 + L 2 ^ 1 ^ 1 ] 



2M W cos/3 



\- [Lf- 1 (m 2 e . tan/3 — M 2 ^ sin 2/3) — L^~ l m ex (jjl + tan/M e J] 



Table 7: Higgs couplings to Sleptons. 



0=1,2 



+2Q 2 igig 2 ig 2 ^(^ 1 ,^ 2 ,s,s)] 



«^ + )l 2 



9 2 X% k {N c -l)\s*pl m \Q%\ 2 
2N C M», 



\Q2^\ 2 \Q^ l \ 2 R(d ka ,s) 



+2QJi Q 22 Qi\ Q12 S(dki,dk2,s,s) 



|M(^4-^°)| 2 = 



g*\>%k{N c -l)\ 
N c M'z cos 2 6»„ 



,a=l,2 

+^?^ (2M 2 (fit - Ml x Ml) + fi 2 s) + (2M 2 (fit - M?M 2 ) + t 2 s) 



+ E ^5 Zjfl g ( g - M D R (^ * 2 ) (Mimi - fi) + fi(s - mi j) 

a=l,2 

+ E ^^ Ks - MljR(u ia , s 2 ) (M 2 (2M 2 a - t) + t(5 - M? J) 



a=l,2 



2|Q|j| S 



(2M 2 (M? i -fi)(M? i -£)-SM) 
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g 2 X'% k (N c -l)\ 
N C M 2 Z cos 2 W 



E i^r'n^T^L^Ka, s 2 ) + 2^i- i ^" 1 ^i l ^ l s 2 ^(d i i, d j2 , § 2 , s 2 ) 



a=l,2 



IQ 



^P^(2Ml{ui-M\Ml) 



+ u 2 s I + 



i^rr 1 ! 2 ^ 



£2 



* {2M 2 z {ut - M\ M 2 Z ) + £ 2 s) 



a=l,2 



■)2i — 1^2*— 1 7 al 7 
«2a V 2 1 ^rf. ^Ufl 



a=l,2 



2|0i*- 1 | a Z. 



§(* - MjjR(d ja , s 2 ) (m 2 (2M? x - «) + - M 2 ^)) 
S(5 - M}jR(d ja , s 2 ) [M 2 {2M 2 jx - t) + £(5 - Mjj) 
^^(2M 2 (M? -m)(M? -£)-stV)l (71) 



where s and p cm are as before. 



C.2.2 Resonant Squark followed by R p Decay 

There are two processes in which a resonant squark is produced via the J3 term in the 
superpotential and these squarks then decay to Standard Model particles: 

1. djd k — s- did m , 

2. uidj — > uid m . 

The matrix elements are given by 

\M.(djd k — > did m )\ 2 = ^Inc' E E ^ / '*ife^ // «'m^ // ™ifc^ // ^"il < 52L| 2 | ( 52/3| 2 (72) 

Q,/3=l,2 i,n=l,3 

5(M ia ,M„^,s,s)s(s-m^ -mjjj 

|M(lijdj — > M;rf m )| 2 = ^jy^ E E ^ / '*ife^ // 'mfc^ // *in^ // ^nl < 52a 1 1 2 1 *^?2^3 ^ (73) 

a,/3=l,2 i,n=l,3 

S(cL, d n/J , s, s)s (s - m 2 ; - m 2 ^) 

C.2.3 Resonant Squark Production followed by Higgs Decay 

There are a number of processes which occur via the production of a resonant squark 
which subsequently decays to either a neutral or charged Higgs. Again we only con- 
sider those processes for which a resonance is possible, i.e. we neglect the processes 
djdk — ► h* R ho/Ho/Ao and u;dj — > d* R ho/Ho/Ao for the first two generations and the 
processes djdk — ► t* 2 ho/Ho/Ao and u ; dj — > b* 2 ho/H /A for the third generation, where 
we consider left/right mixing as these process cannot occur via a resonant diagram: 

1. djdk — > dfH- 

2. djdk — > u^ho/Ho/Ao, 
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4. Ui dj — ► dJih /H /A . 



The matrix elements for these processes are given below. Due to our notation for the 
squark mixing matrices in the case of no left /right mixing the right squark is denoted as 
the second mass eigenstate. 



\M(djd k -> d*H~ 



g 2 (Nc-i)l 

4Af c 



(3 "i2«i/3 

ut - M? M 2 ) + -"^y^ 1 ! 2 ( _ M 2 M 2 



a=l,2 



d i0 H ~ 



dip H ~ 



\M(uidj -> u* k gH 



g 2 (N c -iy. 

4N C 



E ^j^Qfa 1 ? 11 ^^ \ 2 sR(d ka , s) + 2\%Ql k 1 - 1 Ql k 2 - 1 m^- ki m k ^js(d kl , d k 



a=l,2 



iMidjdk^u^H^' 2 



4A"L J |Df| 2 |Q^- 1 | 2 



ui-M 2 - M 2 H+ 



g 2 X% k {N c -l)\ 
4N C 



E IWI^fiaW*-.*) +2£&C&Hl a ^H^iS(ii il ,u a ,S,S) 

a=l,2 

\Q 2 A\ 2 \d\? 



\M( Ui d k ^ d*^ 2 



g2\>% k (N c -l)\ 
4iV c 



Q = l,2 



jldji dj2dj\ 



The coefficients for the various processes can be found in Tables |6] and 
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